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Abstract 


The debate between Anton Tedesko and Charles S. Whitney which 
occurred from the 1930's through the 1950's typifies the confusion among 
designers in the United States regarding thin shell concrete roof design. Each 
man thought his method was correct and designed structures constructed in 
America during the first half of the twentieth century. By taking a closer look at 
their debate, we can gain some insight into their methods of design. To resolve 
the conflict, we then apply modern methods of analysis to analyze a hangar 
model Whitney had presented in his articles. A full span analysis is performed 
using the finite element computer program P-FRAME. In addition, we address 
concerns which were not incorporated into the original analysis. We employ the 
methods of Milo S. Ketchum and Robert S. Rowe to compute deflection 
moments for the structure. In addition, we use Ketchum and Rowe's work as 
background for developing the Initial Deflection Method of computing buckling 
safety factors. To validate the procedure, we compute buckling safety factors 


for a variety of structures and compare them to classical formulations. 
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Chapter One 


Introduction 


Over the last century, designers such as Eugene Freyssinet, Robert 
Maillart, Pierre Luigi Nervi, Felix Candela and Heinz Isler have brought the art 
of reinforced concrete design to its mature state. 

Isler's thin shell concrete roofs cover many European structures. From 
tennis courts to gas stations, his shells provide practical, yet interesting, 
solutions to everyday roofing problems. 

In the United States, however, reinforced concrete design has not 
advanced as it has abroad. A reflection of this lack of progress can be seen in 
the content of basic design texts. During the 1950's, in the well known text book 
by George Winter and Arthur H. Nilson, Design of Concrete Structures, an 
entire chapter was dedicated to arch and shell design. In the 1991 edition, the 
words arch and shell do not even appear in the index.! Why has this type of 
design disappeared from our basic text books? 

A possible explanation is that a general confusion exists in America 
regarding shell behavior and because of this, key safety questions still remain 
unanswered. 

This confusion is clearly demonstrated in the debate between Charles S. 
Whitney (1902-1961) and Anton Tedesko (b. 1902) which took place in the 
1940's and 1950's concerning concrete barrel shell roof design. The design is 
a thin concrete barrel shaped shell with arch stiffeners spaced along the length. 
The debate focused on how to position the shell in relation to the arches. 

Tedesko's opinion was that the shell should be positioned at the rib 


extremity. Whitney, on the other hand, believed that the shell should be located 
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at the mid-height of the rib. They debated in many engineering publications, but 
without resolution. 

Tedesko was an Austrian born engineer who had studied civil 
engineering at the Technological Institute in Vienna in the early 1920's, and had 
learned thin-shell concrete roof design while working at the firm of Dyckerhoff 
and Widmann in Weisbaden. |n 1932, because of prior work experience in the 


U.S., he was sent to work in America when his firm decided to expand its 
operations.^ Once there, he gained an affiliation with the Roberts and Schaefer 


Company in Chicago, and stayed on with them to do extensive thin-shell design 
work during the 1930's and 40's. Because of his effort in this capacity, he 
introduced thin-shell concrete roof structures to the United States.? Examples 
of his work are the Ice Hockey Arena in Hershey, Pennsylvania which opened in 
1936 and the U.S. Navy Hangars at North Island in San Diego, California 
designed in 1941. 

Whitney was an American born designer who wrote two major articles on 
arch design in the 1920's. In the 1940's, he developed a new desing method 


for barrel shell roofs that he called "... a novel feature. . . . which has important 
advantages." 4 |n his method, he placed a great emphasis on volume change 


moments, which are a function of the cross sectional moment of inertia. The 
moments of inertia would be less if the shell was located in the middle of the rib, 
thus, this was the better design. 

Whitney also designed structures which were constructed in the United 
otates. An example of his work is the Field House at Syracuse University built 
in the 1950's. 


We can see, then, that two very different methods of design existed 


‘er = 





simultaneously. The debate between Whitney and Tedesko which started in the 
literature over 40 years ago, has not been resolved. In this thesis, we will first 
examine the different design methods and then utilize modern engineering 
tools to clarify the debate. 

We look at Whitney's method by examining calculations which were 
prepared for a hangar model and presented in articles published during the 
1940's and 1950's. Tedesko's rebuttal to one of the articles is also scrutinized 
to examine his ideas on the subject. 

In our modern analysis, we use the finite element method to analyze 
Whitney's hangar model. Additionally, we employ the methods of Robert S. 
Howe and Milo S. Ketchum to calculate stress amplification due to deflections in 
arches. Using these as a starting point, we develop a method of predicting 


buckling loads for arched structures. 





Chapter Two 
Whitney’s “New Idea” 


Whitney began writing about concrete arch design as early as 1925 with 
his article, “Design of Symmetrical Concrete Arches” published in the American 
Society of Civil Engineers (ASCE) Transactions. Between 1932 and 1940 he 
was the Chairman of the American Concrete Institute (ACI) Committee 312 
which was attempting to establish standards for reinforced concrete arch 
design. His ideas appear in this committee’s reports published in 1932 and 
1940. 

It wasn't until after this that he began writing about his method of 
concrete barrel shell roof design. He presented his ideas in three articles 
published between 1943 and 1955. 

The structure Whitney analyzes in developing his claims is a 220 foot 
clear span aircraft hangar model. The hangar roof is a parabolic barrel shell 
comprised of a 4 inch thick reinforced concrete slab with ribs spaced 20 feet 
center to center. The rib cross section varies from 18 x 32 inches at the crown to 
18 x 40 inches at each springing. The height of the rib center line above the 
supports at mid-span is 27.5 feet, and the roof is supported at each springing by 
concrete A-frames. Whitney assumes that the A-frames act to fix the ends of the 
barrel shell by restricting translational and rotational motion. Figure 1 is a 


longitudinal and transverse section of Whitney's model. This structure was first 
presented in his 1944 article "Aircraft Hangars of Reinforced Concrete".? 


The article that we will focus on was published in the ACI Journal in June 





1950 under the title, "Cost of Long Span Concrete Roof Shells". In this article, 
Whitney explains the advantages of his design method: 


"An important feature of this type of construction is the placing of the shell 
near the neutral axis of the ribs so that the ribs project about half above 
and half below the shell. The principal effects of this arrangement 
structurally are the elimination of edge stresses in the shell due to rib 
flexure and the reduction of the stiffness of the combined rib and shell 


with a corresponding reduction in volume change moments.” € 


He also develops a chart which shows how his shell positioning reduces 


required rib size and thus, construction cost. 


acd e Rıbs 20 ft on center 
Y 


18" y ο crown 


1 


Ribs 2Oft on center 
ALS e 


Ar IB" x 4o" at Springing 
Γιός 


N \ 
| TEN 


| 
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Figure 1. Longitudinal and transverse sections of the 220 ft span hangar model 


His chart, shown in Figure 2, presents three different crown cross 
sections, each based on a 20 foot spacing of the arch ribs. The first cross 
section has the shell located at the mid-height of the rib, the second has the 


shell positioned at the top of a similar rib, and the third has the shell at the top of 
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an enlarged rib. The larger rib in the third cross section is necessary, according 
to Whitney, to provide the equivalent strength of the first cross section. 

To understand Whitney's method of design, we will examine the numbers 
in his chart. Since Whitney did not publish detailed calculations, we must use 
information from his other articles and reports to estimate his results. 

As his first table entry, Whitney gives the moment of inertia at the crown. 
For the first section, it is computed using the 20 foot width and including 
reinforcing steel at the top and bottom of the rib. In the next two sections, he 
only uses a 14 foot width in his calculations. He explains that with the slab at 
the top of the rib, only 7096 of the shell is effective. Whitney does not provide 


any background for this assumption in any of his published material, but we can 


verify it with a formula published in a 1990 textbook on concrete shell design.’ 


The effective overhang of one side of the shell, b,, is: 


1 
be = 0.76 (rhh 


L2 
rT = = , 
where: 8d with: L= Span 
d = Rise 


h = Shell thickness 


Using Whitney’s data for the hangar model: 

be = 6.5 ft 
This gives a total overhang of 13 ft. When the rib width is considered, the total 
effective width becomes 14.5 ft, which compares very well with Whitney’s 14 ft 
assumption. 


Using Whitney’s data at the crown for the hangar model, our computed 


a Y 





moments of inertia for the three cross sections are: 


Section 1: 58,000 in* 
Section 2: 117,200 in^ 
Section 3: 154,800 in^ 


These compare very favorably to Whitney's numbers. Detailed calculations 


n 
7. 

"HN : 

of crown section 58.000 in* 116,300 in- 154,300 in. 

ο... 

live load 1,710,000 in]b.| 1,710,000 in]b.| — 1,710,000 in.lb. 


EXE p ame 
volume change 1,026,000 in.Ib 2,042,000 in.lb. 2 |a 000 in.Ib. 
ῬΜ 3 752,000 in.Ib. όν. 


447,120 Ib. 450,170 Ib. 508,300 Ib. 


Figure 2. Whitney's Data Table 


appear in Appendix A. 











Arch Cross 
Sections 










Total moment 









Maximum horiz- 
ontal thrust 


The second entry is live load bending moment. Whitney uses 30 psf as 
the live load for all three cross sections. For a 20 foot width, this results in a 


distributed load of: 


30 psf x 20ft = 600 2 








In his 1925 article, Whitney derives formulas to compute moments in concrete 


arches for different loads with varying cross sections. From Figure 50 in the 
1925 article, the maximum positive live load moment at the crown is:® 


Mı = Kip L? 


where: p = Uniformly distributed live load 
L = Span length in feet 
K, - Factor from Figure 50 using entering 


arguments N and m 


Ὃν ΡΕ quarterpoint rise - midspan rise 
with: r midspan rise 





I 
m = c with: c = crown values 
I, cos 8, 
s = springing values 


Using Whitney's data for the hangar model: 
m = 0.59 and N =0.25 
With these as entering arguments for table 50: 
K, - 0.0049 
The maximum positive live load moment at the crown is: 


M; - 142,300 ft-Ib. 2 1,710,000 in-Ib 


Whitney presents this value in the table for all three cases. Even though the 
cross sectional length for the second and third cases Is less, he is assuming the 
effective cross section carries the full live load. 

To determine load positioning, Whitney uses influence lines. Figures 35- 
39 in his 1925 article give influence lines for values of "m" ranging from 0.15 to 
0.40 and "N" values ranging from 0.15 to 0.25. From these figures, we can 
extrapolate the load positioning necessary to produce maximum positive 


moment at the crown for m = 0.59 and N = 0.25, which is shown in figure 3. 
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The third value Whitney lists is the volume change moment. The values 


are based on "rib shortening and a temperature drop of 40 degrees F including 


the effect of shrinkage.”? 





| ( 
----------- 0.2} |. — Io L— — Cs — 


Figure 3. Load distribution to produce maximum positive crown moments 


Rib shortening results from compressive axial forces in the structure, and 
the 40 °F temperature drop accounts for the worst case combination of 
temperature change and concrete shrinkage. Obviously, shrinkage produces 
an outward thrust with corresponding negative moments at the supports. With a 
temperature drop, these two effects would add together and make sense when 
considered with the 30 psf live load Whitney is using, which is probably a snow 
load. 

In the "Aircraft Hangars of Reinforced Concrete" article, Whitney provides 
insight to his choice of a 40 °F temperature drop. First of all, he states that 
because plastic flow reduces temperature and shrinkage effects, 60% of the 
maximum temperature range for a geographical area should be used for 
concrete arch calculations. In the article, he presents a table of temperature 


ranges for various locations in the United States. We will choose a value from 





the table for a location in the south, since this is what Whitney implies he used. 
The temperature range of 95 °F for New Orleans, Louisiana is chosen. The 
design range is therefore: 


95 °F x 60% = 57°F 


Assuming that it drops from the mean, our change in temperature will be half of 
the design range, or 27 °F. Whitney suggests that the stress caused by 
shrinkage in concrete can also be represented by a drop in temperature. Using 
his shrinkage value of 15°F results in a total temperature drop 42 °F. 19 
The volume change moments are computed by first calculating the 

thrusts due to rib shortening and temperature changes. The thrusts are 
multiplied by a function of the rise to compute moments. Using Whitney’s 1925 
article, Micalos derives formulas for hingeless (fixed) arches. If we assume a 


secant variation in the cross section, the formula for computing thrust due to rib 


shortening is: 


He. = 45 H le 
RS N 
where: I. = Crown moment of inertia 


h = Mid span rise 
Am = Mean rib area 
H = Dead and live load thrust 


Using Whitney’s data, the thrust due to rib shortening is: 
Hrs = 4,180 lb 


The thrust due to temperature changes is computed from:1? 


I 
Hr = ® aTE= 
ΠΡ 
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where: a = Coefficient of thermal expansion 
T = Temperature change 


Using Whitney’s data and his recommended coefficient of thermal expansion of 
5.5 x10°® in/in/°F, the thrust due to temperature change is:! ? 

H; = 5270 lb 
The crown moment due to the total volume change thrust is: 


M= h (Hrs +H; ) 


all 
3 
With Whitney’s data, the crown moment due to volume changes is: 


M - 1,039,000 in-Ib 


This value compares well with Whitney's value of 1,026,000 in-Ib. Whitney does 
not explain his assuming a secant variation in cross section for this calculation. 

The table entries for the other two sections vary directly with the moments 
of inertia. Thus, the volume change moment at the crown is nearly doubled 
when the shell is shifted from the mid-height to the extremity, and it is increased 
even more for the larger rib. 

The next table entry, the total moment, is simply the addition of the two 
previously calculated values. Since the live load moment is the same for all 
three cross sections, the total moment varies only with the change in volume 
change moments. 

Finally, Whitney computes the maximum horizontal thrust caused by the 


dead and live loads. Whitney's formula for computing thrust due to dead load 


[6:14 
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We b (g - 1) 
rk? 


Hy = 
where: g = The ratio of springing to crown weight 


w. = Crown weight in Ib/ft 
L, = One-half the span length 


k = cosh”! (g) 


Using 150 Ib/ft? as the weight of concrete, the dead load thrust is: 
Ha = 337,600 1b 


The live load thrust for a load that produces the maximum positive moment at 


the crown is:!? 


Cip L? 
H,- B 


where: Cj = A constant computed from Figure 


5lusing m and N 
r - Midspan rise 
p = Distributed live load 


For Whitney’s data: 
H, = 62,800 lb 
The dead and live load total thrust is: 


H = 400,400 lb 


This value is 46,720 lb lower than the table value of 447,120 Ib. Whitney must be 
considering the full span live load in his calculations. To compute the full span 


live load thrust, we will use the previously defined formula with a C; value for 


both maximum positive and negative crown moment live loads:!9 





Πι = 132,000 Ib 


The total thrust would then be: 


H = 469,600 Ib 


This is 22,500 Ib higher than the table value. Even if we consider the negative 
thrust from the volume change calculations, the computed total thrust would still 
be 13,000 Ib higher than the table value. Whitney does not give any explanation 
for this difference. 

In summary, through his table, Whitney shows the importance of volume 
change moments in concrete barrel shell roof design. According to his 
conclusions, placing the shell at the mid-height of the rib cuts the volume 
change moment in half and is the most efficient design. He also claims that with 
the shell located at the extremity, the rib must be increased by 50% to give the 
equivalent strength of a cross section with the shell at mid-height. He does not, 


however, present any calculations to support this claim. 





Chapter Three 


Tedesko's Design ldeas 


Whitney's ideas did not go unchallenged. The discussions of his papers 
raised serious questions as to the validity of his claims. In presenting the 
alternate viewpoint, we will focus on two papers in particular. The first is the 
discussion of Whitney's 1950 paper "Cost of Long-Span Concrete Roof Shells", 
and the second is the discussion of the 1940 report of ACI Committee 312, 
"Plain and Reinforced Arches" . Both discussions were written by Structural 


Engineers from the Roberts and Schaefer Company of New York City, who were 


under the direction of Anton Tedesko at the time.! / Therefore, we will consider 
the alternate viewpoint as Tedesko's. 

Tedesko does not agree with Whitney's claim regarding shell position. 
Tedesko believes, instead, that the shell should be positioned at the rib 
extremity. He develops an alternate chart which shows that the rib width can be 
decreased if the shell is moved to the top or bottom. He supports his arguement 
by looking at the stress distribution and by computing buckling safety factors. 

We will examine Tedesko's chart, Figure 4, to gain more insight into his 
argument. Tedesko uses four crown cross-sections which he calls cases one 
through four. The first two cross sections are the same as presented in 
Whitney's table. The third has the shell at the bottom of the 18 inch wide rib, 
and the fourth has the reduced rib with the shell at the top. 

The first five items are the same ones listed by Whitney. Tedesko carries 
the calculations a bit farther, however, by computing stress distributions, 


tension force taken by the reinforcing steel and a buckling safety factor. 
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Arch Cross 
Sections 


Moment of inertia 

of crown section, 58,000 116,300 116,300 69,700 

in. 

Moment due to E 

live load, in. -1b 1,710,000 1,710,000 1,710,000 1,710,000 
1,026,000 2,042,000 2,042,000 1,230,000 


Concrete stress 
at top of arch 


Concrete stress 
at bottom or arch, 


Total tension 
force in concrete 
to be taken by 


Relative buckling 





Figure 4. Tedesko's data chart 


The moments of inertia for the first two crown sections are the same as 


Whitney's. Tedesko apparently agrees with Whitney's use of the reduced 
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effective width when the shell is moved to the rib extremities. The moment of 
inertia for the third cross section is the same as the second since the two are 
mirror images. Their only differences will be in the section modulus for the top 
and bottom fibers, but this will only affect the stress distribution. For the fourth 
cross section, Tedesko computes the moment of inertia using the full 20 foot 
width. Since the shell is positioned at the top of the rib, however, the reduced 
effective width should be 14 feet. Using the reduced effective width, we 
calculate the moment of inertia to be 65,300 in^; 4,400in* lower than Tedesko's 
table value. 

The moment due to live load is the same as Whitney's and the same for 
all four cases. Tedesko is not challenging Whitney's use of the full width loading 
on the reduced effective width, the amount of load used, or his load position. 

The moments due to volume changes vary directly with the crown 
moment of inertia as in Whitney's table, and the values for the first three cross 
sections are the same as before. Tedesko is not questioning the method of 
computation. The moment for the fourth case is about 6% too high due to 
Tedesko's full width moment of inertia. Using Whitney's volume change 
moment as a starting point, we calculate the moment for the fourth cross section 
to be 1,147,000 in-lb. This is 83,000 in-lb lower than Tedesko's value. 

As in Whitney’s table, the total moments are computed directly from the 
live load and volume change moments. The total for case four using an 
effective width of 14 ft is 2,857,000 in-Ib. 

Tedesko's table values for horizontal thrust for the first three cases are 
the same as Whitney's. He apparently concurs with the calculations from 
Whitney's 1925 article. The thrust for the fourth cross section is computed using 


Whitney's formulas as well. We concur with this value as it is not affected by the 
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change in moment of inertia. 

The next data entry in Tedesko's chart is the stress distribution across the 
crown cross section. This is computed using standard formulas for axial and 
bending stress and the appropriate section modulus. Tedesko assumes an 
uncracked cross section, therefore, an elastic analysis is implied. Since the 
horizontal thrusts are all compressive, and the bending moments at the crown 


are all positive, general formulas for computing stresses at the extreme fibers 


are: 
f = > P = M fo == P + M 
top: A St bottom: A Sp 
where: P = Horizontal Thrust 


A= Cross sectional Area 
M= Total Bending Moment 
S= Section Modulus (I /y) 


Using Whitney’s data for the cross section in case one: 


f, = -1060 psi and fp = + 449 psi 


Tedesko's values are 36 psi lower for both the top and bottom fibers for case 
one. To arrive at the stresses listed in the table for the first cross section, 
Tedesko is either using a total thrust value of 499,800 lb instead of the listed 
value of 447,120 lbs or using a reduced cross sectional area. He does not 
explain the change. The stress calculations for the other three cases also differ 
from what would be expected from Whitney's data, but they do show the trend 
Tedesko wants to demonstrate. In case four, the stresses are almost equally 
distributed and even with the reduced rib, the stresses are still well within the 
strength of concrete. 

Next, Tedesko lists the total tension force to be taken by the reinforcing 
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steel. This is calculated assuming a linear stress distribution and assuming the 
concrete does not carry tension. The first step in computing this value is to 
determine the distance from the bottom fiber to the neutral axis. This can be 


done using a ratio of the table stresses and the rib depth: 


= ae 413 psi M 
Yb = (413 + 1096) psi ρω) ΠΡ 


From this, the total tensile force is: 


T = 5 (413psi)(18 in )(8.75 in) = 32,500 Ib 


This agrees exactly with the table value. With this value for tension, and 
assuming a yield stress in the steel of 50 ksi, the required reinforcing steel area 


would be less than one square inch. The ACI code for minimum reinforcing 


steel area would supersede, thus requiring a reinforcing steel area of:18 


Pmin= 0.005 = as 

solving for A, using the rib cross section yields: 

As = 0.005 (576 in2} = 2.88 in? 
This explains the table reference to minimum required reinforcing steel and 
agrees with Tedesko’s recommended steel area. The tension force and 
required steel areas are calculated similarly for the other cases. For case three, 
the required tension force is zero since the entire cross section is in 
compression at the crown. 


Tedesko computes a safety factor against buckling as the last table value 


for each cross section. The computations are based on Dischinger's formula:!? 
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V. = 3321 El 
H a? 


H = Horizontal Thrust 

E = Effective Modulus of elasticity 
a 

I 


where: 


- One-half the Span Length 
- Moment of Inertia at the Crown 


Using the data for the first cross section with Tedesko's thrust value of 499,820 Ib 
yields a buckling safety factor of: 
MUSS 

This value compares very well with the table value of 8.7. The safety factors for 
the other cross sections are computed using this formula and show how the they 
are increased when the shell is moved to the rib extremity. Even with the 
reduced rib in case four, the safety factor against buckling is larger than when 
the shell is positioned at mid-height. 

Using the same hangar model and data as Whitney does, Tedesko has 
reached the opposite conclusion. In his table, with the shell moved to the rib 


extremity, the rib size can be reduced. Which conclusion is correct? 
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Chapter Four 


Resolving the Conflict 


Our examination thus far reveals questions which must be addressed to 
resolve the shell positioning conflict. In addition to these, other items which 
affect the issue are mentioned in the Whitney and Tedesko articles, but are not 
incorporated into their tables. 

First of all, Tedesko states that additional analysis is required: 


*... the writers do not believe it justified to assume that an investigation 
of the crown of the arch alone can determine the most economic cross 
section. Not only do the maximum moments vary in sign and magnitude 
along the arch axis, but also the relative importance of the volume 
change moments varies. In the lower quarters of the arch the volume 
change moments are only a small percentage of the total design 


moments."20 


Although this is a very serious discrepancy, he does not make the full span 
analysis he claims is necessary. 

Secondly, both sides mention arch deformation effects. And although 
they both imply that the deflections are easily approximated, neither side 
presents any relative data. Whitney even stresses the importance of 


investigating the moments caused by deformations, especially at the crown 
section where "the greatest increase in stress due to deflection occurs."*! 


Tedesko gives some justification for not investigating the additional moments: 


*. . . deformation moments are of important influence only for arches of 
small buckling safety and for arches which do not follow the pressure line 


for dead load" ?? 


He obviously does not consider the deflection moments to be significant in this 


design. 
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Also, the two factions address using a reduced Young's Modulus to 
account for concrete creep when computing deflections. Whitney states that the 


Young's Modulus value for concrete should be reduced by two-thirds to three- 
quarters to calculate deflections under permanent load.^? Tedesko suggests 


using a value of 2,000,000 psi for E. to account for creep.*^ 

To resolve the question of shell position, we will use an approach that 
incorporates analysis methods not available to Whitney and Tedesko and 
addresses the additional points mentioned above. Using the four step process 
outlined below, we will create a table for each of our three cross sections which 
we can use to make comparisons between the two design methods. The 


resulting tables are attached as appendix B. 
REVISED FOUR STEP METHOD OF ANALYSIS 


1 Computing For Moments and Deflection 

We make a full span analysis of the barrel shell roof section utilizing the 
Finite Element computer Program P-FRAME. From the finite element analysis, 
we are able to determine dead load moments, positive and negative moments 
due to different live load distributions, volume change moments due to 
temperature change, axial thrusts at each section and deflections due to the 
loadings. We model three different cross sections and create a table for each. 

We model the arches using the 20 foot width and 4 inch shell thickness 
Whitney specified. The first cross section has the shell at the mid-height of the 
rib and will be referred to as the "Whitney Arch". The second is 14 feet long, has 


the shell located at the lower rib extremity and will be referred to as the 
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“Tedesko Arch". Both of these arches have ribs that vary from 18 x 32 at the 
crown to 18 x 40 at the springing. The third section is 14 foot wide, has a 
smaller rib and has the shell located at the lower extremity. It will be referred to 
as the “Reduced Tedesko Arch". We chose to position the shell at the bottom of 
the rib for the Reduced Tedesko Arch since this will result in a compressive 
Stress distribution through more of the arch span. Since we are designing in 
concrete, this is an important consideration. The rib for this section varies from 
9 x 32 at the crown to 9 x 40 at the springing. The three cross sections are 
shown in Figure 7. 

We used 29 nodes to model each arch, 23 of which are spaced 
horizontally from zero to 220 feet at equal 10 foot intervals. The other six nodes 
are placed to allow for the live loads necessary to make a full span analysis. To 
ensure symmetry, we placed three nodes on each side of the mid-span. Vertical 
positions for the nodes were computed using the equation for a parabola. A 
one-line diagram of the model is shown in Figure 5. 

The moments of inertia at the crown for the first two arches are taken 
directly from Whitney's table. We use the corrected moment of inertia from the 
Tedesxo table for the third. Moment of inertia at the springing is computed using 
Whitney’s dimensions. These calculations appear in Appendix A. To represent 
the varying cross section, we use a linear interpolation between the crown and 


springing, adjusting the value every 10 feet. The Young’s Modulus for all cases 
is 4x 100 psi, the same as used by Whitney and Tedesno, and the coefficient of 
thermal expansion we chose is 5.5 x 107 infin/F, the value that Whitney 


recommends. The end restraints for the arches are modeled as fixed against 


both rotation and translation. For all three cases, we chose a linear elastic 
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analysis. 


To compute dead loads, we input the normal density of concrete, 150 


Ib/ft?, and P-FRAME computes the weights based on cross sectional areas and 
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lengths between the nodes. For the Tedesko arches, an additional externally 
applied dead weight was added to account for the reduced effective width. This 
was modeled as a uniformly distributed horizontal load. The values for this load 


were computed as: 
η] η 2 lb 
(20 ft - 14 ft) ( in) (| 50 ©) = 300 
The live loads are modeled as uniformly distributed and externally 

applied using Whitney’s magnitude of 30 psf. For the full span analysis, we had 
to load the arch with several different live load distributions to produce 
maximum positive and negative moments at the points we wanted to 
investigate. Although data is available through P-FRAME for every nodal point, 
we concentrated our live load analysis on three significant points; the springing, 


the quarter point, and the crown. Each arch is loaded using known distributions 


to produce maximum positive and negative moments at the points of interest.?> 
Load distributions used are shown in Figure 6. 

To model the volume change moments, we applied a uniform 
temperature change of -40 °F along the full span of each arch. The rib 
shortening contribution to the volume change moments is computed directly by 
P-FRAME. The program determines the axial deformations due to the applied 


dead and live loads. 


2 mputing Deflection Moments 
To compute the additional moments due to deflections,we applied the 


theories of Robert S. Rowe and Milo S. Ketchum to the output data from P- 
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FRAME. Both of these procedures are based on series approximations of the 


deflections. To account for the affect of creep, we used a Young's modulus of 


2,000,000 psi in the computer analysis. 
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Figure 6. Live load distributions used for Arch analyses 
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Ketchum's procedure, published in the American Society of Civil 


Engineers (ASCE) Transactions, provides a method for computing final 


deflections as a function of initial deflection, moment and axial force.?9 The 


derivation is based on a beam which is loaded both axially and laterally as 


shown in Figure 7. In the figure, M; is the moment due to q, the lateral load, w; 
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is the deflection caused by the lateral load, w, is the deflection due to the axial 
load and w, is the total deflection. Relating the deflections: 

Wo = Wi + Wa 
Assuming the elastic deflection curves for both loads are similar to their bending 
moment diagrams, a relationship is established between the ratios of the 


deflections and moments at the midpoint of the beam: 








Wa _ Pw 
Wi ML 
Solving for w,: 
τ. απ τ 
Since w, =w,-Ww;: 
Wo Wi + Wi Ed 


Regrouping: 


Solving for w,: 


Thus, the final deflection, w,, can be computed if the initial deflection, the axial 
force and the moment are known. This formula is adapted for our use in 


computing deflection moments for the tables in appendix B as follows: 
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The deflection moment can be computed if the initial deflection, axial force and 
moment at a section are known. The initial data (P, w;, Mj) is available from our 


computer analysis. 


ον -— —— 4744 


Figure 7. Axially and Laterally loaded beam used in the Ketchum derivation. 


Robert S. Rowe's procedure is derrived using a beam loaded both axially 


and laterally as shown in Figure 8. Rowe's method expanded on Ketchum's 
work by applying it to arches as well as beams.?/ Rowe's procedure is once 


again based on the idea of a series of moments. From Figure 8, the deflection, 


Ab, due to the arbitrary lateral load can be expressed as: 


so =o E 


where “n” is a bending moment diagram shape factor. 


27 











The moment caused by this deflection and the given axial load is: 


Mai = P, Ap 





Figure 8. Loaded Beam used in the derivation of Robert S. Rowe's method 


This moment causes an additional deflection of: 





2 
MA Ls EI 


P, ^b L°} 
E I 


At = Mm 


Which, in turn, causes an additional moment: 


2 
M= R A=: P, ny Pade?) 


Assuming the elastic curve maintains the same shape so the n,’s are similar, the 


total moment equation becomes: 
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Multiplying through by [ 24 where A- n, 


Mi(1-A) + P,A s (1- A2) 


M = 
1-A 


With A « 1, A2 goes to zero, therefore: 


Mp[1-A) + PLA 


M = 
J-A 


Assuming the axial load and lateral load bending moment diagrams are similar, 
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Therefore: El 


The total bending moment, including the deformation effects, can be computed 
in terms of the bending moment at the section, the axial load at the section and 
a bending moment diagram shape factor. This is very similar to Ketchum's final 
result. 

Rowe applies this to curved beams. He presents a chart that relates the 
displacement ratios in straight and curved beams to the h/L value in the curved 
beam. From his chart, we see that for arches with rise to span ratios of less 


than 0.15, the deflection in an arch is less than 2% different from that in a 


Straight beam. 


mputing Cros ional Stre 
Stress distributions are computed using the total moments and axial 
thrusts at each section and the standard P/A and My/I stress formulas. Results 


from these calculations appear in each of the tables. 


4 Computing Bucklin fety F r 
A revised method of computing buckling safety factors, based on the 
theories of Rowe and Ketchum, is presented in the following chapter. Data for 


each arch is presented at the end of this chapter. 
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EXPLANATION OF TABULAR DATA 

General 

We show a full span analysis in the three tables in appendix B. Data is 
computed for the crown, the left and right springing and the left and right quarter 
points. The data is displayed in columns from left to right along the arch length. 
An explanation of each line, along with a sample calculation for the "Whitney 
Arch" follows. The P-FRAME output file for the Whitney Arch is attached as 
Appendix C. 


M Bine A 
Moments of inertia for the first two tables are computed based on 
Whitney’s data for the crown and springing rib sizes. The quarter point value is 
linearly interpolated. Data for the third table is calculated based on Tedesko's 

reduced rib width, using Whitney's variation in rib height . 


For Whitney's Arch, the moment of inertia at the crown was previously 


computed as 58,000 in^. For the given rib dimensions of 18 x 40 at the 
Springing, the moment of inertia is 110,000 in^. Using a linear interpolation 


between the two to compute the quarter point value yields 84,000 in. 


Momen nd live | Line (B 

These values are taken from P-FRAME output for the cross sections with 
positive moments acting clockwise at the left and counter-clockwise at the right 
hand end of a segment. The dead load is calculated based on a linear 


interpolation of cross section variation every 10 feet. Dead load for the full 20 
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foot width is applied for all three arches. Load positions for moments due to live 
loads are shown for the four cases in Figure 8. P-FRAME data is converted from 
ft-kips to in-Ib for easy comparison with Whitney and Tedesko data and is 


rounded to four places. For the Whitney Arch at the crown (node 15): 


Moment due to dead load = -30.20 ft-kips = -363,000 in-Ib 


Moment due to live load (1), 
(Maximum positive crown Moment) = +136.4 ft-kips = +1,636,000 in-Ib 


Moment due to live load (2), 
(Maximum negative crown Moment) = -133.7 ft-kips = -1,604,000 in-Ib 


Moment due to live load (3), 

(Maximum positive moment at the 

Left Quarter Point or Maximum 

negative moment at the Left 

Springing) = -64.48 ft-kips 


-774,000 in-Ib 


Moment due to live load (4) 
(Maximum negative moment at the 
Left Quarter Point or Maximum 
positive moment at the Left 


Springing) = +67.10 ft-kips = +805,000 in-Ib 


Initial Displacements, Line (C) 

These values are also taken from P-FRAME output and represent 
displacement of the nodal points from their initial positions due to the indicated 
loads. The displacements shown are in inches and do not include deformation 
moment effects. The dead load displacements are computed using E = 


2,000,000 psi to account for creep. For Whitney's Arch at the crown: 


Initial Displacement due to Dead Load -0.345 inches 
Initial Displacement due to Live Load case (1) = -0.541 inches 
Initial Displacement due to Live Load case (2) = «0.419 inches 
Initial Displacement due to Live Load case (3) = +0.180 inches 
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Initial Displacement due to Live Load case (4) = -0.302 inches 


MEI Line (D: 
Once again, the Values are taken directly from P-FRAME output. The 
thrusts are normal to the indicated cross section and are listed in kips. For 


Whitney's Arch at the crown: 


Axial thrust for Dead Load = 351.6 kips 
Axial thrust for Live load Case (1) = 60.6 kips 
Axial thrust for Live load Case (2) = 70.6 kips 
Axial thrust for Live load Case (3) = 41.1 kips 
Axial thrust for Live load Case (4) = 90.1 kips 


This would give a total thrust of 412,200 lb for dead plus live load for maximum 
positive moment at the crown, and a thrust of 482,800 Ib for dead plus full span 
live load. These values are 3% higher than the thrusts of 400,400 lb and 
469,600 lb computed using Whitney’s formulas. 


Temperature Change Moments, Line (E 

These moments are computed based on the 40°F temperature drop used 
by Whitney applied over the entire arch. The nb shortening contribution, 
computed by Whitney for his volume change moments, is not included here 
since it is computed as part of the Load Moment in line (B). The P-FRAME 
output is once again adjusted from ft-kips to in-Ib for comparison with Whitney 


and Tedesko table values. For the crown section of Whitney's Arch: 


Temperature Change Moment = 59.87 ft-kips = 718,000 in-lb 


If our computed values of dead load, live load for maximum moment at 
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the crown, and temperature change moments are added together, we can make 
a comparison with the total moment value listed in Whitney's table. The P- 
FRAME total moment would be 1,991,000 in-lb. This is 27% lower than 
Whitney's table value of 2,736,000 in-Ib. 


M | ο... 

As the h/L value for our arch is 0.125, according to Rowe, we can 
compute the deflection moments using Ketchum's method. We use P-FRAME 
output to compute the axial thrust, initial deflection and moment at the section for 
each live load condition. The dead and temperature change loads will not 
create deflection moments since they are uniformly distributed across the entire 
arch span. The arch axis will not deflect from the funnicular line under these two 
loads. 

For Whitney's Arch at the crown, the deflection moment at the crown due 


to live load case (1) is computed using: 


M = Pw 
εν 
ML 
From P-RFRAME output 
for Live load case (1) P= 60,600 lbs (compression) 
w; = -0.541 in 


M. - 1,636,000 in-Ib 


The deflection moment at the crown for live load case (1) is: 
M =+ 33,000 in-lbs 
The moment sign is determined by the orientation of the deflection and axial 


force. In this case, the arch is deflecting downward, and the axial force is 
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compressive, thus, a positive moment results. 


Worst Case Moments (G) 

The worst case moments are totaled from Lines B, E and F for dead and 
live loads. The live load case which produces the largest appropriate moment 
is used in each computation. A circled number appears next to the worst case 
moment value to indicate which live load was used. For the crown section of 


Whitney’s Arch: 


Worst case total Moment = +2,024,000 in-Ib 


The worst case total moment for the crown results from a live load for maximum 


positive moment. 


ion Stress, Line (H 
Stresses are calculated based on worst case moments and associated 
axial loads. Standard stress formulas are used with the computed cross 
sectional areas and section moduli for the appropriate point in the arch. For the 
crown section of Whitney’s Arch: 


Axial Load for Worst case Total Moment = 412,200 lb (compression) 


Top and bottom fiber stresses: 


f, = - 840 psi and f, = +277 psi 


Bucklin fety F r 
A safety factor against buckling is computed based on the Initial 


Deflection Method for arches developed in Chapter 5. First, models of the 
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Whitney Arch, the Tedesko Arch and the Reduced Tedesko Arch are given a 
parabolic imperfection. We determine the critical load using the Initial 
Deflection Method, and compute a buckling safety factor by comparing the 
critical load to the dead load plus the 600 Ib/ft live load used by Whitney. The 
finite element analysis program P-FRAME is used to compute the initial 
deflections. 

For each arch, the initial imperfection is a parabola with a midspan rise of 
0.625 ft. We compute the revised nodal coordinates using the transformation 
described in Chapter 5. The revised coordinates are input to the P-FRAME 
program, and we apply a uniformly distributed horizontal load across the entire 
span. To compute the critical load, we must eliminate the rib shortening 
contribution. A purely parabolic model of the arch is loaded with the same 
uniform load, and the resulting deflections are subtracted from those computed 
using the offset model. The final deflections are due to bending moment only. 
The critical load occurs when the deflection is equal to the offset. 

For Whitney’s Arch, at a uniform load of 23.5 kips/ft, the deflections for the 


offset model at the three-quarter point (node 23) are: 
0,73 = -0.224 in and 0,53 7 -2.49 in 
The deflections for the parabolic model are: 
0,73 = -2.76 in and 0,23 = -9.56 in 
This results in a total deflection of: 
6 = -7.51 in = -0.625 ft 


To compute the buckling safety factor, we assume a uniformly distributed dead 
load equal to the average cross sectional value. The buckling safety factor for 


Whitney’s arch is: 
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Safety Factor = 23.5 /2.2 = 10.7 
The P-FRAME input and putput used in this calculation are attached as 
appendix D. Using the Initial Deflection Method for the other two arches yields 
critical loads of: 


Tedesko Arch 47.5 kips/ft 


Reduced Tedesko Arch = — 27.2 kips/ft 


The resulting Buckling Safety Factors are: 


Tedesko Arch 25.0 


Reduced Tedesko Arch 17.4 


These values are higher than the Dischinger values listed in Tedesko's table, 
but show the same trend. Even with a reduced rib, the Tedesko cross section 


has a greater safety against buckling than Whitney's. 
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Chapter Five 


The Initial Deflection Method 
For Computing Critical Buckling Loads 


Using the ideas of Rowe and Ketchum, we will develop a relationship 
between successive deflections and buckling. The general formula is 
developed using an axially loaded column. 

An axially loaded column which is perfectly straight theoretically will not 
buckle, regardless of the applied load. It will only deform along its axis in 


accordance with the well known formula: 


Y 
D 


Òa = 


> 


E 


where: P - Axially applied load 
L = Column length 
A = Cross sectional area 
E - Young's modulus 


Given some type of initial imperfection, however, the column will buckle under 


sufficient load. For the axially loaded column shown in Figure 11 that has an 
initial parabolic offset with a value of 6, at mid-height, the initial moment at the 
mid-span is: 


Mo = P òs 


This moment , in turn, causes an additional deflection, ὃι: 





5 . MyI2 _ (PS )L? 
i 12 EN DIE 
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ΕΠ: 
, B= == - ; 
Letting I2 EI, the additional moment due to the deflection 8, is: 


M¡= Pô = PB 


0.625 in. 





Figure 9. Axially loaded column with initial parabolic offset 


This moment will create an additional deflection , 5,: 


M, L2 
ô= ἜΤ᾽ δοβ' 


The total moment is: 
Mr = Mo+M)+Mo+... = Pd9+P8,+Pd9+...= P dg πο 


Thus, if B is greater than one, the series diverges, the moments will grow without 


bound, and the column will buckle. The critical load occurs when B = 1. When 


B = 1, the original offset, 59, and the initial deflection, 5,, will be related as: 
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δι Ξ B 09 = do 


We have defined the critical point in terms of deflection. When the initial 
deflection is equal to the original offset, the column is at the critical point. If the 
initial deflection is greater than the original offset, the column will buckle. This 
criteria is easily applied to output from finite element computer programs. To 
check buckling, one only needs to compare the computer generated deflections 
to chosen input imperfections. 

As acheck, we will compare the results from the “Initial Deflection 


Method” with several well known buckling formulations. 


Euler Column Buckling 


We will check the initial deflection buckling criteria against the classic buckling 


problem presented by Euler. His solution for the critical load of a column 


hinged at both ends is: °9 


Por = nel 
L2 
where: P = Column axial load 


L = Length between the supports 
E = Young’s Modulus 
I = Minimum moment of inertia of the cross section 


The physical model used is a 100 inch long steel beam with a 12 square inch 
cross section which is 6” x 2 “. The beam is hinged at both ends. 


Since the chosen cross section results in a minimum moment of inertia of 


4 in^, the critical load using Euler's formula is: 
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Per = 118 kips 


To test the Initial Deflection Method, the beam is modeled on the finite 
element computer program SAP-90 using frame elements with 11 equally 
spaced vertical nodes. The horizontal offset for each node is calculated using a 
parabolic equation with the chosen maximum offset as 0.625 inches at the mid- 
span. The column is modeled with a pin at the bottom and a roller at the top. 
Loads are applied in the negative vertical direction at the top node. The critical 
load is one that produces a 0.625 inch deflection at the middle node. To 
determine the critical point, we start at the theoretical critical load and perform 
iterations until we read a deflection of 0.625 inches on the output. 

Using the Initial Deflection Method with SAP-90, the critical load is : 

Per = 116 kips 


This is 1.7% lower than Euler’s theoretical value. 


Plate Buckling 


Next, we check the initial deflection buckling criteria against the theory for a 
simply supported plate. The physical structure we model is a 100 inch by 100 
inch steel plate which is one inch thick. 

The general formula for critical load per length for a simply supported 


plate uniformly compressed in one direction is:3° 





(Nx) cr xs τ (B+ af 


where: a= Horizontal plate dimension 
b = Vertical plate dimension 
N, = Load per unit length along the horizontal 
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DC Eh 
and: 12 (1 - v2) 


where: E = Young’s modulus 
h = Plate thickness 


v = Poisson’s ratio 
For a square plate, a = b, and this becomes: 


2 
(Να) - AD 
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For the chosen plate, with E 2 30 x 106 psi and v - 0.3, the flexural rigidity is: 
D = 2,747 in-kips 


From the Timoshenko formula, the critical distributed load is: 


kips 


(Nie = 1084 TE 


The structure is modeled using the finite element computer program SAP-90 
with 100 plate elements. The computer model is composed of 121 nodes with 
11 nodes spaced equally along the horizontal, and 11 vertical nodes equally 
spaced at each of these. Each vertical line of nodes has a parabolic offset with 
a maximum of 0.625 inches at the middle node. The plate is simply supported 
on all sides with a pin along the bottom edge and rollers along the other three. 
We model the uniform load across the top of the plate as a pressure load. 

As with the axially loaded column, the critical loading occurs when the 
mid-span horizontal deflection equals the original offset. Using this criteria, we 
load the plate at the theoretical critical load, and adjust as necessary until we 
see a mid-span deflection of 0.625 inches in the computer output. 


Using the Initial Deflection Method, the critical load is: 
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(NJa= 911 Sins 
in 


This result is 15.6 % lower than the theoretical buckling load. This result is 


suprising in light of the success with the column data. 
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Arch Buckling 


In a two hinged parabolic arch, uniformly distributed loads are carried 
axially to the supports. Bending moments are zero throughout, and the arch 
simply "squats" due to rib shortening. By introducing an initial imperfection as 
we did with the axially loaded column, however, bending moments are created 


(Figure 12 refers). The moments can be expressed as: 
Mo- P do 


where: P = Axial load at the section 
9, = initial offset from the funicular line 
The greatest moments will occur at the points where the largest offset is located. 
These bending moments will, in turn, cause additional deflections. These 


deflections are expressed as they were for the axially loaded column as: 


5 - ΜοΙ2 (Ρδ͵}12 
= JEI 12E1 








Using the same derivation as in the case of the axially loaded column, we can 
define the critical point of arch buckling in terms of deflection. When the initial 
deflection is equal to the original offset, the arch is critically loaded. If the initial 
deflection is greater than the original offset, the arch will buckle. To compute 
buckling loads using the Initial Deflection Method, we need to compare 
computer generated deflections to our chosen offsets. 

To validate the Initial Deflection Method for arch buckling, we will 
compare it with Timoshenko's theoretical results. The model we use is a 2 
hinged steel arch with a constant 6" x 2" cross section. The arch spans 100 


inches horizontally and has a rise of 10 inches. 
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Timoshenko's general formula to compute critical loads for a uniformly 


loaded parabolic arch with a constant cross section is:3! 


Ga = Aa El 
13 
where: E - Young's modulus 


] = Moment of inertia 
L = Horizontal span length 


Ay = A factor depending on the height to span 
ratio and the number of hinges 
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Figure 10. Bending Moments in Arches caused by Innial Offset 





From Timoshenko's Table 7-5 with h/L 2 0.1, 4, 2 28.5. The critical load with E 


- 30 x 100 is: 
Gcr = 3.42 kips per inch of horizontal span 
The Initial Deflection Method is tested using the finite element computer 
program P-FRAME. We model the arch using 21 nodes equally spaced along 
the horizontal and compute initial vertical coordinates for these nodes using a 
parabolic equation with a rise of 10 inches at the middle node. The arch is free 
to rotate and restricted from translating at both ends. 


We will apply a parabolic offset to each half span as shown in Figure 11. 





Figure 11. Parabolic Offset used in the Computer Model 


The maximum offset will be 0.625 inches at each quarter point. Since we want 
the nodes to be offset from the initial parabolic curve, each node will have to be 


adjusted both vertically and horizontally as shown for node 5 in Figure 12. The 
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change in coordinates will be a function of the parabolic offset, ὃς, and the angle 
Os as: 


ôy 5 = ὃς (cos Qs) and Ó, 5 = ὃς (sin Qs) 





Figure 12. Blown up view of the coordinate transformation at node 5 


The arch is loaded using uniformly distnbuted honzontal loads. To 
eliminate the affect of nb shortening from the computer output, a purely 
parabolic arch model without initial offsets is loaded with a uniformly distnbuted 
Critical load. The nb shortening deflections are subtracted from the deflections 
computed using the offset model. The resulting deflections are due to bending 
moment only. Once again, the theoretical buckling load is chosen as a Starting 
point and iterations are performed until we see an adjusted deflection of 0.625 
in either half of the arch. Using the Initial Deflection Method, the buckling load 


for the arch model is: 
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Ρο = 2.92 EDS 
in 


This is 14.6% lower than Timoshenko’s theoretical value of 3.42 kips/in. 
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Chapter Six 


Conclusions 


From the debate between Charles S. Whitney and Anton Tedesko we 
can draw several conclusions. First of all, Whitney's 1925 article was a good 
guideline for arch design. His formulas and charts give information which is 
backed up by modern methods of analysis. Using only graphic statics and the 
Calculus, he computed formulas and charts to determine thrusts and moments 
for significant loads and load positions at several points in the arch. For 
maximum positive moment at the crown, his distribution is less than 196 different 
from accepted design guidelines used from the 1950's to today. His adjustment 
factors allow quick computations for a very diverse range of arch designs. 
Although his conclusions for thin shell design are vehemently challenged by 
Tedesko, his method of calculating arch thrusts, live load moments and volume 
change moments are not. 

Secondly, Tedesko's claim that a full span analysis is necessary for arch 
design is valid. By looking at our full span analysis, we can see that stress 
distributions at the springing, quarter point and crown must all be investigated. 
The crown analysis that Whitney and Tedesko both present is not sufficient to 
design a barrel shell roof. 

Next, both men were correct in their claim that the deflection moments 
are not significant for Whitney's hangar model at the assumed live loads. We 
can see, however, that the deflection moments and corresponding stresses 
must be addressed in barrel shell roof design. As the factor of safety against 


buckling is reduced, the deflection moments become more significant. The 
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methods of Rowe and Ketchum are a reliable way to determine deflection 
moments from data available from most finite element programs. 

Additionally, Whitney's emphasis on volume change moments is valid. 
His claim that the shell should be located at the mid-height of the rib to reduce 
volume change moments is also valid. There is an irrefuteable relationship 
between the moment of inertia and the volume change moment. When we look 
at stress distributions over the full span and take all loads into account, 
however, we see that the decreased moment of inertia has a drawback. 
Heducing the moment of inertia increases the tensile stresses near the 
supports. The tensile stress at the springing for Whitney's Arch is double that of 
Tedesko's Arch. 

Also, there is a significant advantage in placing the shell at the bottom of 
the rib as opposed to the top. Because concrete is weak in tension and strong 
in compression, we want to reduce tensile stresses as much as possible. 
Placing the shell at the bottom of the rib does this. We can see from the stress 
distributions for Tedesko's Arch and Tedesko's Reduced Arch that tension only 
exists in upper part of the cross section for the first and last quarter of the arch. 

The Initial Deflection Method is a tool which can be used with existing 
finite element computer programs to compute buckling safety factors for a 
variety of structures. Although longhand computations of initial offset values are 
tedious, this could easily be written into a finite element computer program. The 
critical loads computed with the Initial Deflection Method are conservative for 
plates, but could provide a ready check for buckling capacity for arches during 


an iterative design process. 
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APPENDICIES 


APPENDIX A: CALCULATIONS FOR MOMENTS OF INERTIA 
APPENDIX B: DATA TABLES FOR THE FOUR STEP METHOD OF ANALYSIS 
APPENDIX C: SAMPLE COMPUTER DATA FILE FOR WHITNEY'S ARCH 


APPENDIX D: SAMPLE COMPUTER DATA FILE FOR INITIAL DEFLECTION 
METHOD CALCULATIONS 





Moment of Inertia Calculations 


ion 1: 18” thick rib with shell | mid-heigh 

As = As’ = 3.0 inches, n = 7 

Io A Yt Ayt d to na 
Rib 49,200 576 16 9220 0 
Left Shell 592 444 16 7100 0 
Right Shell 592 444 16 7100 0 
Top Steel D Zl 2:5 35 1575 
Bottom Steel 2 21 29.5 620 1955 

y 2 16.0 in I = 58,000 in* 
ection 8” thick ith shell located at top (Flange 70% effective 

As = As’ = 3.0 inches, n = 7 

I A γι Ay, d to na 
Rib 56,800 618 16 9890 6.9 
and Steel 
Left Shell 400 300 2 600 τι 
Right Shell 400 300 2 600 Fol 

y =9.1 in I = 117,200 in4 


26 


Ad? 


3830 
3830 


Ad? 
29.400 


15,100 
15,100 





ion 3: 27" thick rib with shell | 


As = As’ = 4.5 inches, n = 7 


I A 
Rib 73,700 864 
Left Shell 376 282 
Right Shell 376 282 
Top Steel 3 32 
Bottom Steel 3 32 
y= 10.7 in 


ion 2: 9" thick rib with shell | 


As = As’ = 1.5 inches, n = 7 


ls A 
Rib 24,580 288 
Left Shell 450 318 
Right Shell 450 318 
Top Steel 0 1.5 
Bottom Steel 0 1.5 
y = 25.6 in 


at top (Flanae 7095 effective 
Yt Ay, d to na 
16 13,800 5.3 
2 564 8.7 
2 564 8.7 
255 80 8.2 
29.5 944 18.8 
I = 154,800 in4 
70% effective 
Yt Ay, d to na 
16 4608 9.6 
30 9540 4.4 
30 9540 4.4 
A 3 21.6 
30 15 4.4 
I = 65,300 in* 
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24,300 
24 500 
21,300 

2,150 
11,300 


Ad? 


26,740 
6160 
6160 
700 
30 





Whitney Arch 


ett Left Crown Right Right 
Springing Qtr Pt QtrPt Springing 
A) Moment of Inertia 110,000 84,000 58,000 84,000 110,000 
(in4) 
B) Load Moment (x10 in-ib) 
Dead Load -1.936 +0.321 -0.363 . 40.321 -1.936 
Live Load 1 :2.575 -1.535 41.636 -1.535 +2.576 
Live Load 2 -2.812 41.492 -1.604 1.492 -2.812 
Live Load 3 -6.572 +2.996  -0.774  -1.686 44.159 
Live Load 4 6.335 -3.010 0805  À 41.643 -4.396 
C) Initial Displacements (in) 
Dead Load 0 -0.476 -0.345  -0.476 0 
Live Load 1 0 40.196 -0.541 +0.196 0 
Live Load 2 0 -0.260 «40.491 -0.259 0 
Live Load 3 0 +1.155 0.180 -0.917 0 
Live Load 4 0 41.092  -0.302 40.980 0 
D) Axial Thrust (kips) 
Dead Load 397.4 362.3 351.6 362.3 397.4 
Live Load 1 61.8 62.8 60.6 62.8 61.8 
Live Load 2 85.0 725 70.6 72.5 85.0 
Live Load 3 56.8 42.9 41.1 41.2 39.7 
Live Load 4 89.9 92.4 90. 1 94.0 107.1 
E) Temperature Change -1.774 40.096 +0.718 +0.096 -1.774 
Moment (x 10° in-1b) 
F) Deflection Moments 
Live Load 1 0 -0.012 40.033 -0.012 0 
Live Load 2 0 40.019 -0.030 40.019 0 
Live Load 3 0 40.050 -0.007 -0.039 0 
Live Load 4 0 -0.104 10.028. «0.098 0 
G) Worst Case Total 
Moment (x106 in-1b) -10.38 (3) +3.433(3) +2.024(1) +2.158(4) -8.106(4) 
H)Section Stress (psi) 
Top +1593 -996 -840 -7595 +1168 
Bottom -2144 +475 +2// +169 -1780 
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Tedesko Arch 


Eeft Left Crown Right Right 
Springing Qtr Pt Qtr Pt Springing 
A) Moment of Inertia 218,400 167,400 116,300 167,400 218,400 
(in4) 
B) Load Moment (x10 in.1b) 
Dead Load -2.354 0.268 40.002 «0.268 -2.354 
Live Load 1 +2.459 -1.527 41.687 -1.527 -2.459 
Live Load 2 -2.941 +1.497 -1.555 +1.497 +2.941 
Live Load 3 -6.643 +2.972 -0.744 -1.685 +4.080 
Live Load 4 +6.160 -3.002 10.876  Á 41.655 -4.562 
C) Initial Displacements (in) 
Dead Load 0 -0.474 -0.663  -0.474 0 
Live Load 1 0 40.078 -0.309 +0.078 0 
Live Load 2 0 -0.155 +0165. 5155 0 
Live Load 3 0 -0.595 +0.064 +0.446 0 
Live Load 4 0 40.518 -0.208  -0.524 0 
D) Axial Thrust (kips) 
Dead Load 400.5 365.5 354.5 365.5 400.5 
Live Load 1 61.3 62.3 60.1 62.3 61.3 
Live Load 2 84.5 71.9 70.0 71.9 84.5 
Live Load 3 56.6 42.6 40.7 40.9 39.4 
Live Load 4 89.3 91.7 89.4 93.3 106.4 
E) Temperature Change -3.506 +0.193 +1.424 +0.193 -3.506 
Moment (x106 in-1b) 
F) Deflection Moments 
Live Load 1 0 -0.005 «0.019 -0.005 0 
Live Load 2 0 40.011 -0.012 +0.011 0 
Live Load 3 0 +0.026 -0.003 -0.018 0 
Live Load 4 0 -0.048 +0.019 «0.050 0 
G) Worst Case Total 
Moment (x106 in-1b) -12.50 (3) +3.481(3) +3.132(1) +2.166(4) -9.939(4) 
H)Section Stress (psi) 
Top +1261 -844 -958 -684 +897 
Bottom -1028 -96 -95 -219 -923 
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Reduced Tedesko Arch 


ΡΕ Left Crown Right Right 
Springing Qtr Pt Qtr Pt Springing 
A) Moment of Inertia 123,800 94,600 £65,300 94,600 123,800 
(in4) 
B) Load Moment (x10 in-ib) 
Dead Load -1.506 40.135 40.001 40.135 -1.505 
Live Load 1 +2.509 -1.532 +1663 -1.532 -2.509 
Live Load 2 -2.891 41.496 -1.573  Á 41.496 42.891 
Live Load 3 -6.617 42.970 -0.755  -1.684 44.113 
Live Load 4 416.235 -3.005 40.846 41.649 -4.495 
C) Initial Displacements (in) 
Dead Load 0 -0.507 -0.715 -0.507 0 
Live Load 1 0 40.153  -0.520 40.153 0 
Live Load 2 0 -0.256 10.325 -0.256 0 
Live Load 3 0 -1.041 240.132 40.800 0 
Live Load 4 0 40.938  -0.327  -0.903 0 
D) Axial Thrust (kips) 
Dead Load 318.9 291.7 282.9 291.7 318.9 
Live Load 1 61.5 62:5 60.3 62.5 61:5 
Live Load 2 84.7 72.2 70.2 72.1 84.7 
Live Load 3 56.7 42.7 40.9 41.0 39.6 
Live Load 4 89.5 92.0 89.7 93.6 106.7 
E) Temperature Change -1.989 +0.107 +0.805 0.107 -1.989 
Moment (x10 in-1b) 
F) Deflecuon Moments 
Live Load 1 0 -0.010 | 40.032 -0.010 0 
Live Load 2 0 0.019 -0.023 40.019 0 
Live Load 3 0 0.045 -0.005 -0.033 0 
Live Load 4 0 -0.089 «0.030 «0.089 0 
G) Worst Case Total 
Moment (x 10® in-Ib) -10.11 (3) +3.257(3) +2.501(1) +1.980(4) -7.989(4) 
H)Section Stress (psi) 
Top +2189 -1355 -1350 -1001 +1600 
Bottom -1074 -94 -126 -245 -978 
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CRITICAL LOAD ON WHITNEY ARCH 
ΧΧΧ xxx 


| 7 | τι bute dl lead 
m Sloped UDL Proj. UDL Local UDL Local UDL Triangular Thermal 
= K/ft slope K/ft horiz k/ft perp K/ft parll K/ft @ GJ Change (F) 


1 ο 223,3 0 ο ο 0 
2 0 223.9 0 ο 0 ο 
3 0 5-5 0 0 0 0 
ü ο 3069555 ο 0 0 ο 
5 ο 23235 0 ο O 0 
6 ο 22339 O 0 O 0 
7 ο 223,0 0 ο 0 0 
8 ο 5 0 0 O 0 
p ο 5-5 0 ο 0 0 
Ὁ 0 mz. 0 0 0 0 
E 0 m2349 0 O 0 0 
p ο 55 0 0 0 O 
B ο 223.5 0 0 0 ο 
4 ο ο... 5 0 0 0 ο 
I5 O 55 0 ο 0 0 
6 ο = Ce O 0 0 ο OF: 
p 0 E.» 0 0 ο 0 
8 0 E239. 0 0 0 0 
p 0 =23.9 0 0 ο O 
0 ο 2555 0 0 ο 0 
E ο 223,5 0 ο ο ο 
p 0 223,5 0 0 0 0 
B 0 223,9 0 ο 0 0 
4 ο 25.59 0 0 0 0 
Io 0 55 0 ο ο 0 
6 O E24. 0 ο 0 0 
7 0 EO 3. ο O ο ο 
8 0 m2 39 ο ο 0 0 


ped UDL, Projected UDL & Point Loads act in the global coordinate system. 
tal Perpendicular, Local Parallel, Triangular Loads act in 

| local member coordinate system. 

angular Loads are O at the lower joint with the magnitude specified at 

> greater joint. 


ΕΥ 





= | 


Joint Load 
Number Case 


r 1 
2 1 
3 1 
4 1 
> i 
6 1 
7 1 
8 1 
E 1 
10 1 
11 1 
12 1 
13 1 
l4 1 
15 J 
16 al 
17 T 
18 1 
19 1 
20 I 
Pl 1 
22 1 
23 1 
24 i 
25 Í 


WHITNEY'S ARCH 


X-Displ. 
(in) 


. 00000 
.04866 
-01909 
.05002 
„12296 
19793 
. 22448 
. 24696 
126596 
. 24262 
-2 3091 
.20001 
. 18682 
2102 29 
.00000 
210299 
.18682 
«20001 
223001 
. 24262 
226996 
. 24696 
~ 22448 
2197903 


212996 
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Y-Displ. 
(in) 


0.00000 
-.16462 
=. 91029 
ποτ Ἐ 
-1456679 
=. 17921 
"2953 
πο ου ο” 
-3.41666 
-3.93464 
-4.04874 
-4.28527 
-4.36562 
-4.66788 
cu 17595 
-4.66788 
-ᾳ. 36562 
-4.28527 
-4.04874 
-3.93464 
=3.4 1666 
S2 
= 2.485833 
ο. 11εἰ 


56517 


CRITICAL LOAD ON PARABOLIC ARCH 


xxx 


Rotation 
(rad) 


0. 


00000 


.00183 
.00320 
.00413 
.00470 
.00482 
.00506 
.00520 
.00460 
.00372 
.00366 
.00331 
.00322 
.00171 
.00000 
.00171 
.00322 
.00331 
.00366 
.00372 
.00460 
.00520 
.00506 
.00482 


.00470 





CRITICRESEOAD ON PARABOLIC ARCH 


WHITNEY'S ARCH 


ase Results 


Joint Load X-Displ. Y-Displ. Rotation 

lumber Case (ia) SOT, (rad) 
26 1 05002 ὉΠ Τα .00413 
27 1 OSOS ος.» 200/320 
20 D .04866 5 160662 -00 T83 
29 i ο οσο ο 0.00000 0.00000 
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ase Results 


oint 
lumber 


1 


2 


10 


1 


12 


n3 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


Load 
Case 


1 


a 


CRITICAL LOAD ON WHITNEY ARCH 


xxx 


X-Displ. 


T2: 


E 


22 


(in) 

.00000 
233012 
«89038 
.45051 
.87802 
. 12684 
.18403 
„13963 
.14320 
.03244 
.00495 
.95114 
. 93647 
721131 
2199977 
.18664 
. 43349 
. 48345 
251526 
.66847 
292652 


183337 


76211 
63772 


21944 
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Y-Displ. 


0. 


VO 


EO 


300 


(in) 


00000 


. 40387 


ο το 


.48996 


. 39780 


.87432 


89919 


. 74565 


-96761 


«99980 


-16961 


.10760 


31791 


-61165 


.01194 


223630 


.00986 


-29158 


92319 


13806 


90458 


206335 
«96067 
«88741 


«08717 


xxx 


Rotation 
(rad) 


0 


.00000 
. 00688 
00962 
200213 
.00628 
2000795 
.00086 
200373 
.00885 
"01357 
.01479 
„07710 
01787 
201909 
201 9659 
201708 
ED di 
-01129 
.00805 
.00654 
.00018 
200/679 
.00969 
"01237 


.01691 


Ae -. -- —— Ram 





'ase Reau Ite 
Joint Load 
lumber Case 


È 


a 





CRITICAL LOAD 


X-Displ. 
( in) 


pp 592398 
-.86867 
-.23086 


0.00000 


ON WHITNEY ARCH 


Y-Displ. Rotation 
(in) (rad) 
το 652861 -01892 
22261023 -01736 
-.77844 201128 
0.00000 0.00000 











CRITICAL LOAD ON WHITNEY ARCH 


xxx ΧΧΧ 
X - coord. Y - coord. X - Degree Y - Degree 
(feet) (feet) of Freedom of Freedom 
0 0 0 0 
9.919 4.963 1 1 
19.87 9.438 1 1 
29.84 132043 1 1 
39.83 16.92 1 1 
49.84 19.92 1 al 
54.86 21.24 1 1 
59.88 22.42 1 1 
69.92 24.43 1 1 
79.94 25.95 l 1 
82.55 26225 1 1 
87.96 26.8 l JE 
89.98 26.96 1 1 
100 27.48 1 1 
110 24.5 1 1 
120 27.08 1 1 
129.98 26.22 1 1 
131.98 26 1 1 
137.46 ο. οι 1 1 
139.95 24.97 1 d 
149.94 23.31 T 1 
159.92 21022 1 1 
164.88 20.02 1 1 
169.86 187/2 1 1 
179.83 15.82 1 1 
189.89 12.49 1 1 
199.87 8.774 1 1 
209.99 4.583 1 1 
220 0 0 0 


Degree of Freedom: O=restrained 1=free j=coupled to joint 
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